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Mathematical analysis of plasmonic nanoparticles: the scalar case^ 
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Abstract 

Localized surface plasmons are charge density oscillations confined to metallic nanopar¬ 
ticles. Excitation of localized surface plasmons by an electromagnetic field at an incident 
wavelength where resonance occurs results in a strong light scattering and an enhancement 
of the local electromagnetic fields. This paper is devoted to the mathematical modeling of 
plasmonic nanoparticles. Its aim is threefold: (i) to mathematically define the notion of 
plasmonic resonance and to analyze the shift and broadening of the plasmon resonance with 
changes in size and shape of the nanoparticles; (ii) to study the scattering and absorption 
enhancements by plasmon resonant nanoparticles and express them in terms of the polar¬ 
ization tensor of the nanoparticle. Optimal bounds on the enhancement factors are also 
derived; (iii) to show, by analyzing the imaginary part of the Green function, that one can 
achieve super-resolution and super-focusing using plasmonic nanoparticles. For simplicity, 
the Helmholtz equation is used to model electromagnetic wave propagation. 
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1 Introduction 

Plasmon resonant nanoparticles have unique capabilities of enhancing the brightness of light and 
confining strong electromagnetic fields |35|. A thriving interest for optical studies of plasmon 
resonant nanoparticles is due to their recently proposed use as labels in molecular biology [22] . 
New types of cancer diagnostic nanoparticles are constantly being developed. Nanoparticles are 
also being used in thermotherapy as nanometric heat-generators that can be activated remotely 
by external electromagnetic fields m- 

According to the quasi-static approximation for small particles, the surface plasmon res¬ 
onance peak occurs when the particle’s polarizability is maximized. Plasmon resonances in 
nanoparticles can be treated at the quasi-static limit as an eigenvalue problem for the Neumann- 
Poincare integral operator, which leads to direct calculation of resonance values of permittivity 
and optimal design of nanoparticles that resonate at specified frequencies 121 El ED Eo] ED. At 
this limit, they are size-independent. However, as the particle size increases, they are determined 
from scattering and absorption blow up and become size-dependent. This was experimentally 
observed, for instance, in [23l [331136] . 

In |6], we have provided a rigorous mathematical framework for localized surface plasmon 
resonances. We have considered the full Maxwell equations. Using layer potential techniques, we 
have derived the quasi-static limits of the electromagnetic helds in the presence of nanoparticles. 
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We have proved that the quasi-static limits are uniformly valid with respect to the nanoparti¬ 
cle’s bulk electron relaxation rate. We have introduced localized plasmonic resonances as the 
eigenvalues of the Neumann-Poincare operator associated with the nanoparticle. We have de¬ 
scribed a general model for the permittivity and permeability of nanoparticles as functions of the 
frequency and rigorously justified the quasi-static approximation for surface plasmon resonances. 

In this paper, we first prove that, as the particle size increases and crosses its critical value for 
dipolar approximation which is justified in [H], the plasmonic resonances become size-dependent. 
The resonance condition is determined from absorption and scattering blow up and depends on 
the shape, size and electromagnetic parameters of both the nanoparticle and the surrounding 
material. Then, we precisely quantify the scattering absorption enhancements in plasmonic 
nanoparticles. We derive new bounds on the enhancement factors given the volume and electro¬ 
magnetic parameters of the nanoparticles. At the quasi-static limit, we prove that the averages 
over the orientation of scattering and extinction cross-sections of a randomly oriented nanopar¬ 
ticle are given in terms of the imaginary part of the polarization tensor. Moreover, we show 
that the polarization tensor blows up at plasmonic resonances and derive bounds for the ab¬ 
sorption and scattering cross-sections. We also prove the blow-up of the first-order scattering 
coefficients at plasmonic resonances. The concept of scattering coefficients was introduced in 
[S] for scalar wave propagation problems and in [TO] for the full Maxwell equations, rendering 
a powerful and efficient tool for the classification of the nanoparticle shapes. Using such a 
concept, we have explained in |3| the experimental results reported in mi. Finally, we con¬ 
sider the super-resolution phenomenon in plasmonic nanoparticles. Super-resolution is meant 
to cross the barrier of diffraction limits by reducing the focal spot size. This resolution limit 
applies only to light that has propagated for a distance substantially larger than its wavelength 
uniiiTi. Super-focusing is the counterpart of super-resolution. It is a concept for waves to be 
confined to a length scale significantly smaller than the diffraction limit of the focused waves. 
The super-focusing phenomenon is being intensively investigated in the field of nanophotonics as 
a possible technique to focus electromagnetic radiation in a region of order of a few nanometers 
beyond the diffraction limit of light and thereby causing an extraordinary enhancement of the 
electromagnetic fields. In mm, a rigorous mathematical theory is developed to explain the 
super-resolution phenomenon in microstructures with high contrast material around the source 
point. Such microstructures act like arrays of subwavelength sensors. A key ingredient is the 
calculation of the resonances and the Green function in the microstructure. By following the 
methodology developed in mm, we show in this paper that one can achieve super-resolution 
using plasmonic nanoparticles as well. 

The paper is organized as follows. In section we introduce a layer potential formulation 
for plasmonic resonances and derive asymptotic formulas for the plasmonic resonances and 
the near- and far-fields in terms of the size. In section we consider the case of multiple 
plasmonic nanoparticles. Section is devoted to the study of the scattering and absorption 
enhancements. We also clarify the connection between the blow up of the scattering frequencies 
and the plasmonic resonances. The scattering coefficients are simply the Fourier coefficients 
of the scattering amplitude 13 no). In section we investigate the behavior of the scattering 
coefficients at the plasmonic resonances. In section[^we prove that using plasmonic nanoparticles 
one can achieve super resolution imaging. Appendix]^ is devoted to the derivation of asymptotic 
expansions with respect to the frequency of some boundary integral operators associated with 
the Helmholtz equation and a single particle. These results are generalized to the case of multiple 
particles in Appendix [B| In Appendix [C] we provide the technical modifications needed in order 


3 


to study the shift in the plasmon resonance in the two-dimensional case. In Appendix we 
prove useful sum rules for the polarization tensor. 

2 Layer potential formulation for plasmonic resonances 

2.1 Problem formulation and some basic resnlts 

We consider the scattering problem of a time-harmonic wave incident on a plasmonic nanopar¬ 
ticle. For simplicity, we use the Helmholtz equation instead of the full Maxwell equations. The 
homogeneous medium is characterized by electric permittivity Em and magnetic permeability 
/im, while the particle occupying a bounded and simply connected domain L) d (the two- 
dimensional case is treated in Appendix of class for some 0 < a < 1 is characterized 
by electric permittivity Ec and magnetic permeability /ic, both of which may depend on the 
frequency. Assume that ?R.Ec < 0, ^Ec > 0, ?R.fic < 0, Sl/Uc > 0 and define 

km — 


and 

= EmX{^\D) + EcX{D), HD = EmX{^^\D) + £cX(^), 

where x denotes the characteristic function. Let u*(x) = be the incident wave. Here, a; 

is the frequency and d is the unit incidence direction. Throughout this paper, we assume that 
Em and Hm are real and strictly positive and that ^kc < 0 and Qkc > 0. 

Using dimensionless quantities, we assume that the following set of conditions holds. 

Condition 1. We assume that the numbers Em, He cire dimensionless and are of order 

one. We also assume that the particle has size of order one and u is dimensionless and is of 
order o(l). 

It is worth emphasizing that in the original dimensional variables oj refers to the ratio between 
the size of the particle and the wavelength. Moreover, the operating frequency varies in a small 
range and hence, the material parameters Ec and Hc can be assumed independent of the frequency. 
The scattering problem can be modeled by the following Helmholtz equation 




V -—Vu + uF‘EDU = t) in 
HD 

u+ — u- = 0 on dD, 


1 du 
Hm dv 


1 du 
He dv 


= 0 


on dD, 


u® := u — satisfies the Sommerfeld radiation condition. 


( 2 . 1 ) 


Here, d/dn denotes the normal derivative and the Sommerfeld radiation condition can be ex¬ 
pressed in dimension d = 2, 3, as follows: 


du 

d\x 


'ikm.U 


< C\x 


(d+l)/2 


as Ixl —>■ -|-oo for some constant C independent of x. 
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The model problem (2.1) is referred to as the transverse magnetic case. Note that all the 
results of this paper hold true in the transverse electric case where ed and are interchanged. 
Let 


Fi{x) = —u\x) = —e 


ikmd-x 


1 /)?/^ 1 

F2{x) = - = - kme^^-<^-^d-u{x) 

OV firn 

with i^(x) being the outward normal at a: G dD. Let G be the Green function for the Helmholtz 
operator A + satisfying the Sommerfeld radiation condition. In dimension three, G is given 
by 


G{x,y,k) = -- 


Jk\x-y\ 


47r|x — y\ 

By using the following single-layer potential and Neumann-Poincare integral operator 


•Sd['^]{x)= G{x,y,k)ip{y)da{y), 
JdD 


X G 


{jc%rmx) = 


dG{x,y,k) 


fan d^{x) 

we can represent the solution u in the following form 


'4’iy)dcr{y), X G dF>, 


u(x) = 


u^ + S 


n 


X G 


X G D, 


( 2 . 2 ) 


where ^ H 2 [dD) satisfy the following system of integral equations on dD [7]: 

=Fi, 

I + =F2, 


(2.3) 


where Id denotes the identity operator. 

We are interested in the scattering in the quasi-static regime, i.e., for w <C 1. Note that for 
(jj small enough, is invertible [7]. We have cj) = [S^[ijj] — Fi), whereas the following 

equation holds for ip 

Ad{uj)['iP] = f, 


where 


1 /I . . 1 /I 


Ad{uj) = —{-Id+{}C%-)*) + -{^Id-{F>S)*){S'pp)-^S'pp, 


fj>m 2 


flc 2 


1 /I 


/ = F2 + -{-Id-iJC>p,r)iS'pp)-\F,]. 

Me ^ 


(2.4) 

(2.5) 

( 2 . 6 ) 


It is clear that 


yiD(o) = Ad.o = T(1 h + +1 (1h-= (^ + T),d_ (i _ ijx:;,. ( 2 . 7 ) 

fAm ^ f^c ^ -^Mc Me t^m 
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where the notation 1C\) = (/C^)* is used for simplicity. 

We are interested in finding We first recall some basic facts about the Neumann- 

Poincare operator ]C*q 

Lemma 2.1. (i) The following Calderon identity holds: ICd<Sd = SdK,*j-,; 

(a) The operator }C*jj is self-adjoint in the Hilbert space H~^{dD) equipped with the following 
inner product 

{u,v)'H* = -{u,Sd[v])_ii ( 2 . 8 ) 

with (•, •)_! 1 being the duality pairing between H~^{dD) and H^idD), which is equivalent 
2 ’ 2 

to the original one; 

(in) Let 'H*{dD) be the space H ~2 (dD) with the new inner product. Let (Xj, ipj), j = 0 , 1 , 2 ,... 
be the eigenvalue and normalized eigenfunction pair of ICfy in 'H*{dD), then Xj G (— 5 , 5 ] 
and Xj —)• 0 as j ^ 00 ; 

(iv) The following trace formula holds: for any ip G 'H*{dD), 


(—-/d + lCf))[ip] 


dSppip] 

dn 


(v) The following representation formula holds: for any ip £ H 


j=0 


It is clear that the following result holds. 

Lemma 2.2. Let ^-[{dD) be the space H^{dD) equipped with the following equivalent inner 
product 

{u,v)'H = {{-Sd)~^[u],v)_ii. ( 2 . 9 ) 

Then, Sp is an isometry between 'H*{dD) and ^{{dD). 

We now present other useful observations and basic results. The following holds. 

Lemma 2.3. (i) We have {—^Id JCf))Sf)^[x{dD)] = 0 with x{dD) being the characteristic 

function of dD. 

(a) Let Ao = 5 , then the corresponding eigenspace has dimension one and is spanned by the 
function (po = cSfy^[x{dD)] for some constant c such that Hf/Joll-H* = 1- 

(Hi) Moreover, 'H*{dD) = h ^ C}, where H^idD) is the zero mean subspace 

ofn*{dD) and pj G H-UdD) for j > 1, i.e., {pj,x{dD))_i^i = 0 for j > 1. Here, {pj}j 
is the set of normalized eigenfunctions of Kf). 
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From (2.7), it is easy to see that 


where 


Ti = 



(2.10) 

j=0 


^ 1 2 1 1 ^ 

2 /ic he hm 

(2.11) 


We now derive the asymptotic expansion of the operator A{uj) as cj —>• 0. Using the asymp¬ 
totic expansions in terms of k of the operators S^, and (/C^)* proved in Appendix [a| 

we can obtain the following result. 

Lemma 2.4. As a; —)> 0, the operator Ad{oj) : 'H*{dD) —)■ 'H*{dD) admits the asymptotic 
expansion 

,2 


where 


Ad{^^) = Ada + ^^Ad,2 + 0{uj^), 

A ^ \ir I ^rnhm — £chc J J r'* \ c -1 ( 

Ad, 2 = \£m- £c)I^D,2 + - 


Me 


Id-lCo)S^^SD,2- 


( 2 . 12 ) 


Proof. Recall that 


(2.13) 


Ad(^) = —(h<i + (Kj;")-) + 

Mm ^ l^c ^ 

By a straightforward calculation, it follows that 

{S^) = Id + u{^A^cIJ‘C^D,iSd A y/£mhmS£)^SD,l) + 

{^chc^0,2^0 + V £chc£'mhm^D,lSD,l + ^mhmSj^SD,2) + 0{ljJ^), 

— Id + Smhm y/Acflc^<S£) Sjo i + 

^ S£) 2 “b ^chd^/^chc a/S jo iSjj 'Sr),l) 

+0(a;3). 


where Bd,i and Bd ,2 are defined by (A.5). Using the facts that 


and 


{-Id-JCh)S^^SD,i = 0 


^Id - {JC’hT = {\ld - ICh) - eiCD,2 + o{k^), 


the lemma immediately follows. 


□ 


We regard Ad(w) as a perturbation to the operator A_d,o for small oj. Using standard pertur¬ 
bation theory |34] . we can derive the perturbed eigenvalues and their associated eigenfunctions. 
For simplicity, we consider the case when Xj is a simple eigenvalue of the operator /CJj. 

We let 

Rji = ("^£>,2[<7’j], (2-14) 
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where Ad ,2 is defined by (2.12). 

As oj goes to zero, the perturbed eigenvalue and eigenfunction have the following form: 


where 


Tj{uj) 


= 


= 


= Tj -\- Uj‘^Tj^2 + 0{UO^), 

(2.15) 

= ifj -b uAtpj^2 0{uj^), 

(2.16) 


(2.17) 


(2.18) 


2.2 First-order correction to plasmonic resonances and field behavior at the 
plasmonic resonances 

We first introduce different notions of plasmonic resonance as follows. 

Definition 1. (i) We say that uj is a plasmonic resonance if 

|rj(a;)| <C 1 and is locally minimal for some j. 


(a) We say that uj is a quasi-static plasmonic resonance if \Tj\ <C 1 and is locally minimized 
for some j. Here, tj is defined by (2.11). 

(Hi) We say thatuj is a first-order corrected quasi-static plasmonic resonance if\Tj-\-oj'^Tj^ 2 \ ^ 1 
and is locally minimized for some j. Here, the correction term Tj ^2 is defined by (2.17). 


Note that quasi-static resonances are size independent and is therefore a zero-order approx¬ 
imation of the plasmonic resonance in terms of the particle size while the first-order corrected 
quasi-static plasmonic resonance depends on the size of the nanoparticle (or equivalently on uo 
in view of the non-dimensionalization adopted herein). 

We are interested in solving the equation .4 .d(w)[</>] = / when uo is close to the resonance 
frequencies, i.e., when Tj{u}) is very small for some j’s. In this case, the major part of the solution 
would be the contributions of the excited resonance modes (pj{uj). We introduce the following 
definition. 


Definition 2. We call J C N index set of resonance if tj’s are close to zero when j € J and 
are bounded from below when j G H. More precisely, we choose a threshold number 770 > 0 
independent of uo such that 

\Tj\ >Vo>0 for j £ H. 

Remark 2.1. Note that for j = 0, we have tq = IjfJm, which is of size one by our assumption. 
As a result, throughout this paper, we always exclude 0 from the index set of resonance J. 

From now on, we shall use J as our index set of resonances. For simplicity, we assume 
throughout that the following conditions hold. 

Condition 2. Each eigenvalue Xj for j £ J is a simple eigenvalue of the operator IC)^. 






Condition 3. Let 


(2.19) 


^ _ f^m “T l^c 

2 (ltm /^c) 

We assume that X ^ 0 or equivalently, fic ^ —fJ-m- 

Condition 1^ implies that the set J is finite. 

We define the projection Pj{oj) such that 

In fact, we have 

Pj{u) = Y,P^{u:) = Y,i^. I {^-AD{uj))-^d^, ( 2 . 20 ) 

jeJ jeJ 

where 'jj is a Jordan curve in the complex plane enclosing only the eigenvalue Tj(u}) among all 
the eigenvalues. 

To obtain an explicit representation of Pj{lj), we consider the adjoint operator 
By a similar perturbation argument, we can obtain its perturbed eigenvalue and eigenfunction, 
which have the following form 


Tj{ui) = Tj{uj), (2.21) 

(pj(uj') = ipj-\-uj‘^(pj^2pO{Lo‘^). ( 2 . 22 ) 

Using the eigenfunctions we can show that 

Pj{uj)[x\ = (x, (2.23) 

j&J 

Throughout this paper, for two Banach spaces X and Y, by C{X, Y) we denote the set of 
bounded linear operators from X into Y. 

We are now ready to solve the equation AD{oj)[iX\ = /• First, it is clear that 

^ = + AD{u;)-HPMu:)[f]]. (2.24) 

In 


The following lemma holds. 

Lemma 2.5. The norm ||.Td(w)~^T’j=(w)||£('H*(9 D),'H*(aD)) uniformly bounded in oj. 

Proof. Consider the operator 

AD{io)\jc : Pj.{uj)n*{dD) ^ Pj.{u:)n*{dD). 

For u small enough, we can show that dist(iT(.4,£)(a;)|jc),0) > where a{AD{’pj)\j<^) is the 
discrete spectrum of Ad{<jj)\j‘^. Then, it follows that 

||ylD{w)-‘(Pj.(<.>)/)ll = ||(.4 d(i..)|p,.)“‘(Fj.(w)/)|| < iexp(il)||Pj.H/||, 

% ^0 
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where the notation A < B means that A < CB for some constant C. 
On the other hand, 

j&J j&J 


teJ 


Thus, 


\\PAu;)\\ = \\{Id-Pj{com<{l + 0{u;)), 

from which the desired result follows immediately. 


□ 


Second, we have the following asymptotic expansion of / given by (2.6) with respect to uj. 

Lemma 2.6. Let 

fi = (-[d- t^ix)] + —{]-Id-IC}j)S]^^[d- {x - z)]^ 

V L'm Me ^ / 

and let z be the center of the domain D. In the space P*{dD), as oj goes to zero, we have 

/ = w/i + 0{u?), 

in the sense that, for oj small enough, 

\\f-u^h\\n*<Cu:^ 


for some constant C independent of oj. 

Proof. A direct calculation yields 

/ = F2 + U\ld-{lC)fr){S^^T\Fi] 

t^C ^ 

= — • l/(x)] + O{0j‘^) + 

f^m 

^{il^d - )Ch) {{Sd)-^ + ujBd,i) + O{oj^)) l-P'^-^-^ixidD) + " ^)]) + 0{oj^)] 

pikm.d‘z 1 , .pikmd'Z i 

=- —{^Id-KDS^^xidD)] - ——[-ld-lC*j,)BD,i[x{dD)] - 

ooiyte^H^e"^’-^'^-^ (—[d ■ n{x)] + —{^Id-IC}))Sf)^[d- (x - + 0(a;^) 

\ Me ^ J 

= -uiy/£mtimP^”^‘’‘'^ (—[d-v{x)] + —{^Id - lC}))Sff)\d ■ (x-z)]^ 

V hm Me ^ / 

+0(w^), 

where we have made use of the facts that 

{hd-lC*j,)S-^\x{dD)]=t) 
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and 


Bn,i[x{dD)] = cS^%{dD)] 

for some constant c; see again Appendix [A) □ 

Finally, we are ready to state our main result in this section. 

Theorem 2.1. Under Conditions\^ and[^ the scattered field = u — vh due to a single 
plasmonic particle has the following representation in the quasi-static regime: 


where 




^ - +0{u}), 


jeJ 


Tj{u) 


2 ^ -^- X , - + 0(w) 


J6J 


\-\j+ 0 (^ 2 ) 


with A being given by (2.19). 
Proof. We have 


; {f , 4 I- x-l/D ( \f\ 

tp = 2 ^-- +Ad{u:) iPjfiu})f), 


KJ 




+ 2 ^ - (tZ - + 0 (^ 2 ) 


0 (a;). 


jG.J 2/i, 


We now compute (/i, with /i given in Lemma 2.6 We only need to show that 

(^{\ld-ICD)Sfi^[d-{x-z)]),:pj^^ ={d-v{x),ipfiH*- (2.25) 
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Indeed, we have 


- }C*D)S^^[d • {x - z)], 


n* 


- {S^^id • {x - z)], {hd - }Cn)SD[^j]) , 

- ( 5^1 [d-{x- z)],Sd {\ld - ICh) 

-(d-{x- z), j 


= -(d-{x-z),- 


dSni^j 


IdD 


d[d ■ {x — z)] 
dv 


du 


Sni^Pjjda 


1 1 
■ 2 ’ 2 


- J • (x - z)]SD[^j] - ASD[(fj][d • (x - z)]^dx 


where we have used the fact that is harmonic in D. This proves the desired identity and 

the rest of the theorem follows immediately. □ 


Corollary 2.1. Assume the same conditions as in Theorem 2.1. Under the additional condition 
that 


min |Tj(a;)| S> uj , 
j&J 


(2.26) 


we have 


V’ = 2^-^^-^7^-TT^i-+ 0{uj) 


X-Xj+uj^il-J^) r,-2 

More generally, under the additional condition that 


min To (w) ^ w 
j&J 


m+1 


for some integer m > 2, we have 


i’ = Y. 


ikr, 


^ikmd-z 


{d ■ u{x), (pj)'u*Tj + 0(a;2) 




ieJ 


-1 


+ 0(a;). 


ij,m 


Rescaling back to original dimensional variables, we suppose that the magnetic permeability 
/ic of the nanoparticle is changing with respect to the operating angular frequency u while that 
of the surrounding medium, is independent of w. Then we can write 




(2.27) 


Because of causality, the real and imaginary parts of obey the following Kramer-Kronig 
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relations; 


+ 00 


1 


/(a;) 

n'{uj) 


TT 


-p.v. 


J —oo 
r*+oo 


CU — S 

1 


-fi'{s)ds, 


TT 


-p.v. 


UJ — s 


^ {s)ds, 


(2.28) 


where p.v. stands for the principle value. 

The magnetic permeability can be described by the Drude model; see, for instance, 

|35] . We have 

= Ho{l - F — - 2 I ■ -1 )’ (2-29) 

where r > 0 is the nanoparticle’s bulk electron relaxation rate (r“^ is the damping coefficient), 
T is a filling factor, and ojq is a localized plasmon resonant frequency. When 


(1 - F)(a;2 - urlf - Fuj^ioo'^ - w^) + r"V < 0, 


the real part of //c(<^) is negative. 

We suppose that D = z + 5B. 
that 


The quasi-static plasmonic resonance is defined by w such 

1^ ^J‘ra + _ < 

2(/rm -/^c(w)) ^ 


for some j, where Xj is an eigenvalue of the Neumann-Poincare operator It is clear 

that such definition is independent of the nanoparticle’s size. In view of (2.15), the shifted 
plasmonic resonance is defined by 


argmm 


2 /J"n 




where Tj 2 is given by (2.17) with D replaced by B. 


3 Multiple plasmonic nanoparticles 

3.1 Layer potential formulation in the multi-particle case 

We consider the scattering of an incident time harmonic wave tt* by multiple weakly coupled 
plasmonic nanoparticles in three dimensions. For ease of exposition, we consider the case of L 
particles with an identical shape. We assume that the following condition holds. 

Condition 4. All the identical particles have size of order S which is a small parameter and the 
distances between neighboring ones are of order one. 

We write Di = zi + 6D, I = 1, 2,..., L, where D has size one and is centered at the origin. 
Moreover, we denote Dq = 5D as our reference nanoparticle. Denote by 

L 

D = [JDi, £o= £raX{^\D) + £cX{D), fiD = PmX{^^\D) + ficXiD). 

l=l 
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The scattering problem can be modeled by the following Helmholtz equation: 


V • — Vtt + uP'eDU = 0 in M?\dD, 
IJ-D 

u+ — U- = 0 on dD, 


1 du 
Hm du 


1 du 
/Uc du 


= 0 on dD^ 


u® := u — satisfies the Sommerfeld radiation condition. 


Let 


u\x) = 

Fi,iix) = -u\x)\Q^^ = -e 


ikmd-x I 


\dDU 


^ , du'', , 

Fi,2{x) = 


dD, 


= • u(x)\g^ , 


and define the operator by 


Fdp,Di['^]{x) — 


I a Dp 


dG{x,y,k) 

du{x) 


'4’iy)dcr{y), x G dDi. 


Analogously, we define 

•SDp,D,['^]ix) = [ G{x,y,k)'tp{y)da{y), x G dDi. 


I a Dp 


The solution u of (3.1) can be represented as follows: 

L 


u{x) = < 


a"+ 

i=i 

Esljw. 


x G M.^\D, 
x G D, 


1=1 


where (t)i,il>i G H '^{dDi) satisfy the following system of integral equations 

•^ix.D, i'^pi — Fi,i, 

p^i 

i{hd + (rj,7)-)|*] + Gild- {ic%y)m 

l^m ^ h^c ^ 

f^m ,, 
p^l 


(3.1) 
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and 


Fi^i = -u* on dDi, 

1 

Fi ,2 = - on dDi. 

Hm. du 

3.2 First-order correction to plasmonic resonances and field behavior at plas- 
monic resonances in the multi-particle case 

We consider the scattering in the quasi-static regime, i.e., when the incident wavelength is 
much greater than one. With proper dimensionless analysis, we can assume that a; <C 1. As a 
consequence, is invertible. Note that 


- Fi,,). 


p¥=i 

We obtain the following equation for V’i’s, 

Ad{vj)[iP] = /, 

where 



( Ad^ (w) 

\ 


( ° 

-4i,2(w) 

-4l,L(w)^ 

Ad{w) = 

AD 2 X 



-4,2,1 (ca) 

0 

-42,l(w) 


+ 



0 : 


V 

Adi^ x) 




-4l,l_i(w) 0 / 


V’ = 


/V’A 

■02 

V0 l/ 


, / = 


//i\ 

/2 

V/A 


and 


1 A 


Ap(^) = -{5'‘'-(CD,)*)(5D,)-'5Dr.D. + —KfT.D,. 

fAc ^ fAm 

fi = Fi,2 + U\ld-{lC%r){S^j,^X\Fi,,]. 

Pc ^ 

The following asymptotic expansions hold. 

Lemma 3.1. (i) Regarded as operators from 'H*{dDp) into T-L*{dDi), we have 

AdjX = Aojfl + 

(ii) Regarded as operators from 'H*{dDi) into 'H*{dDj), we have 

AiX) = - {\ld - XXo] (5z,p,o,i + Xa2) + + 0(5"). 

Me ^ l^m 
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Moreover, 


{hd-}Ch,)oS^^oSi^p,0,i = o{5^), 
{hd - JCh,) o O Si,p,0,2 = 0 { 6 ^), 

ICi,p,o,o = 0(<5'). 


Proof. The proof of (i) follows from Lemmas 2.4 and B.3 


We now prove (ii). Recall that 


= lld-)Ch,+0{6^u;% 

= S^l-k,SilSD„iSil + Oi6W), 

^^,Dp ~ ‘^^,p,0,0 + ‘5z^p,o,l + Si,p,0,2 + kmSl,p,i + k‘^Sl,p,2,0 + 0{6^) + 0(w^(5^) 

^d:,a = }Ci,p,o,o + 0{u2^6^). 


Using the identity 

{hd-}Ch,)S^^[x{Di)] = 0, 

we can derive that 

- }C*jj^){S'ffJ ^S’fffjj +—}Ci,p,o,o + 0{6'^uj'^) 
t^c ^ ^ h^m 

= fC-id- Kd,)s^‘s‘7o + + 0{SV) 

h^c ^ l^m 

— — “ ^*Di)S£)^ (Sl,p,o,0 + + Si,p,0,2 + kmSi,p,i + k^Si,p,2,0 + 0{5^)) 

/ic Z 

H- ^i,p,o,o + 0(S‘^co'^) 

l^m 

— — ~ ^dJSd^ {Si,p,o,i + Si,p,o, 2 ) H- K,i,p,o,o + 0{5^uP‘) + 0{6^). 

The rest of the lemma follows from Lemmas IB.31 and IB.61 □ 

Denote by 'H*{dD) = 'H*{dDi) x ... x 'H*{dDL), which is equipped with the inner product 

L 

i' 4 ’, 4 >)n* = 4 >i)H*(dDi)- 

i=i 


With the help of Lemma |3.1[ the following result is obvious. 

Lemma 3.2. Regarded as an operator from %*{dD) into T-L*{dD), we have 

> 4 ( w ) = j^D,o + -^D,! + 0(10^6“^) + O(S^), 
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where 


b 

O 

II 

/ Aoifl 

Ad2 ,0 

\ 

, Ad,i = 

/ 0 

Ad,1,21 

Ad, 1,12 AD,l,i:i 

0 Ad,1,22, 

... \ 

1 

V 

j 


\ Ad,i,li 

Ad,1,LL-1 

0 / 


with 


•^D,lpq 


'2Hm ^Hc 

1 ,l 


1 


~ ~ ^*Dp)<SdI {•^P,Qfl,l + <Sp,q,0,2) + —)<^p,q,0,0- 

f^c ^ l^m 


It is evident that 


where 


00 L 

EE 

3=0 l=l 


— + — - f—- —)A- 

2/im 2/ic /ic Hm 

(fj^l = (fjCl 


(3.2) 


(3.3) 

(3.4) 


with ei being the standard basis of M^. 

We take A{u}) as a perturbation to the operator Aofi for small co and small 5. Using a 
standard perturbation argument, we can derive the perturbed eigenvalues and eigenfunctions. 
For simplicity, we assume that the following conditions hold. 

Condition 5. Each eigenvalue Xj, j £ J, of the operator JCf)^ is simple. Moreover, we have 

< b. 

In what follows, we only use the first order perturbation theory and derive the leading order 
term, i.e., the perturbation due to the term Ad,i- For each I, we define an L x L matrix Ri by 
letting 


h^l,pq — ) 

= (^AD,i[piep],g3ie^^^, 

= {■^D,l,pqVPl]-,^l)y^.>- 

Lemma 3.3. The matrix Ri = {Ri^pq)p^q=i^...^L has the following explieit expression: 


Ri.,pp 

Rl,pq 




15 x°‘y^(pi{x)(pi{y)da{x)da{y) 

\^p ^q\ 


MTTyUc 

0(<53), 



pAq- 



X * y 

- nPi{x)pi{y)da{x)da{y) 

Zp - Zq]-^ 
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Proof. It is clear that Ri^pp = 0. For p ^ q, we have 


Rl,pq — Rl,pq + Rljpq + Rljpq-, 


where 


R 




nil 

,pq 


flc 

1 /A 


jlc 

1 


(^( 2 -Id lCp)^Sj^^Sp^qfl^2{Pl\iPl 


n*idDi) 


n*{dDi) 


^l,pq {l^p,q,0fi[pi]: Pl)qi*(dDi)' 

We first consider Rjpq- By the following identity 

did-)Ch)SDApi]=SD,{\id-)CD,)[pi] = (A, - hifi, 


we obtain 


d^i,pq ~ ^ ^Dp)‘SDl‘^P,qfi,l[Pl]:‘SD, 


L^dDi) 


Using the explicit representation of 5p^q^o,i and the fact that {x{dAj), (pi)L‘2{dDj) = 0 for J A 0) 
we further conclude that 


d^i,pq — 0- 


Similarly, we have 


nil 

^l,pq 




l^c 




Hzp - ^ 

JdDq JdDo 47r|2;p - Zq\^ 4Tr\zp - Zq\^ 


\a\=m=l 


')pi{x)(pi{y)do{x)da{y) 


Airy, 


o' I I 15 x°'yd ^pi{x)^pi{y)da{x)da{y) 

^ \a\ = \l 3 \ = l'>9DoJdDo \Zp Zq\ 


O'l-l) E 


+ v^(Ai - J) E / / I- ^ — ^x'^y'^ipi{x)ipi{y)da{x)da{y). 

2 Jgjoq JqDo \zp-Zq\-^ 


Finally, note that 


where Um 


^p,qfi,oYPl\ — 


rO • u{x) = 


dTrlzp-ZqP dvrlzp-ZqP 

((y Zq')ni, Pi)i^ 2 (^Q£)^y and a = (cii, (I 2 , (I 3 ) • 


3 

^ ^ O-ml'mix), 

m=l 
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By identity (2.25), we have 


^l,pq ^ (^P:'?>0,o[¥^«]) V^O'H*(aDi) 


4vr|zp 

-Zq 

P/^m 


1 


dvr Zp 

-Zq 



1 


47r|zp 

-Zq 



1 


4vr|zp 

-Zq 

P/^m 




H*{dDi) 


L^(aDp) 

X ■ yipi{x)ipi{y)da{x)d(j{y). 


(A, - T.) 


I dDo J a Do 


This completes the proof of the lemma. 


□ 


We now have an explicit formula for the matrix Ri. It is clear that Ri is symmetric, but not 
self-adjoint. For ease f presentation, we assume the following condition. 

Condition 6. Ri has L-distinct eigenvalues. 

We remark that Condition is not essential for our analysis. Without this condition, the 
perturbation argument is still applicable, but the results may be quite complicated. We refer to 
|26j for a complete description of the perturbation theory. 

Let Tj^i and Xj^i = ^Xj^i^iY', I = 1,2, ...,L, be the eigenvalues and normalized 

eigenvectors of the matrix Rj. Here, T denotes the transpose. We remark that each Xj^i may 
be complex valued and may not be orthogonal to other eigenvectors. 

Under perturbation, each Tj is splitted into the following L eigenvalues of ^(w), 

Tj^i{u}) = Tj + Tjj + + O(co^d^). (3.5) 

The associated perturbed eigenfunctions have the following form 

L 

X 0{6^) + 0{uP‘5‘^). (3-6) 

p=i 

We are interested in solving the equation Ad{^)['4’] = / when cj is close to the resonance 
frequencies, i.e., when Tj{u)) are very small for some fs. In this case, the major part of the 
solution would be based on the excited resonance modes ^pjj{(jj). For this purpose, we introduce 
the index set of resonance J as we did in the previous section for a single particle case. 

We define 

13 / \ / \ / j ^ Ji 

= I ^ J., 

In fact, 

Pj(u;) = ^Pjiuj) = E ^ - ^Diuj))-^d(, (3.7) 

j&J j&J 
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where 7 ^ is a Jordan curve in the complex plane enclosing only the eigenvalues for I = 

1,2,... ,L among all the eigenvalues. 

To obtain an explicit representation of Pj{uj), we consider the adjoint operator Ad{oj)*. By 
a similar perturbation argument, we can obtain its perturbed eigenvalue and eigenfunctions. 
Note that the adjoint matrix .RJ = Rj has eigenvalues tJJ and corresponding eigenfunctions 
Xj^i. Then the eigenvalues and eigenfunctions of Ad{oj)* have the following form 

= Tj+Tj;i + 

+ 0(5^) + 0 ( 0 ;^ 

where 

L 

p=i 

with Xj^i^p being a multiple of Xj^i^p. 

We normalize ipj^i in a way such that the following holds 

(,^j,p: ^j,q)'H*idD) — ^pqi 

which is also equivalent to the following condition 


^j,p ^j,q ~ ^pq- 


Then, we can show that the following result holds. 

Lemma 3.4. In the space 'H*{dD), as oj goes to zero, we have 

/ = w/o + 

where fo = (/o,i, • • •, /o,l)^ with 

fo,i = (—d-iy{x) + — • (x - = 0{6^. 

V hm he ^ / 

Proof. We first show that 

\\'^\\'H*(dD(i) = W'^W'U-^dOy W'^W'HidDo) = W'^W'HidD) 


for any homogeneous function u such that u{Sx) = 5^u{x). Indeed, we have r]{u){x) = S^u{x 


Since ||??(^i)ll-H*(aD) = ^ ^ lkllw*(aDo) (see Appendix 

3 

\MnRdDo) = ^'^\\viu)\\n*{dD) 


B), we obtain 


= ji+™| 


u\ 


H*(dDy 


which proves our first claim. The second claim follows in a similar way. Using this result, by a 
similar argument as in the proof of Lemma 2.6 we arrive at the desired asymptotic result. □ 


Denote by Z = [Zi,..., Zi), where Zj = ■ We are ready to present our main 

result in this section. 
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Theorem 3.1. Under Conditions 010 an the scattered field by L plasmonic particles in 
the quasi-static regime has the following representation 


where 


^ ^— Ti I I 

jeJ l=l ^ ' 

,3 


~ S A - A, + (^ - ^)-V,,, + 0(^4) + OiS^u;^) 
Proof. The proof is similar to that of Theorem |2.1[ 

As a consequence, the following result holds. 


□ 


Corollary 3.1. With the same notation as in Theorem 3.1 and under the additional condition 
that 


for some integer p and q, and 


min |t,';( a;)| S> 


Tj,l{uj) = Tj^l^p^q + o(u;'^dP), 


we have 


^ (d ■ Tj)'H*{dDo)^^j,l Tj,l + 0{uj‘^d2) ^ 


jeJ 1=1 


T,i 


j,cp,q 


4 Scattering and absorption enhancements 

In this section we analyze the scattering and absorption enhancements. We prove that, at the 
quasi-static limit, the averages over the orientation of scattering and extinction cross-sections of 


a randomly oriented nanoparticle are given by (4.10) and (4.11), where M given by (4.7) is the 
polarization tensor associated with the nanoparticle D and the magnetic contrast pifioj)/iXrn- In 


view of (4.15), the polarization tensor M blows np at the plasmonic resonances, which yields 


scattering and absorption enhancements. A bound on the extinction cross-section is derived 


in (4.17). As shown in (4.20) and (4.22), it can be sharpened for nanoparticles of elliptical or 


ellipsoidal shapes. 


4.1 Far-field expansion 

For simplicity, we assnme throughout this section that D contains the origin. We first prove the 
following representation for the scattering amplitude. 

Propsition 4.1. Let x G be such that |x| 3> 1/uj. Then, we have 




^ikm\x\ 

47 r|x| 




(4.1) 
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with 


A. 




IdD 


being the scattering amplitude and V' being defined by (2.3). 


Proof. We recall that the scattered field can be represented as follows: 

u^{x) =S^j^[fi]{x) 


1 

dvr 


'dD 


^ikm\x-y\ 

\x - y\ 


'fiiy)dcr{y). 


From 


it follows that 


\x-y\ = \x\(l- + O(r^) ) , 


u^{x) = — 


(4.2) 




which yields the desired result. 


□ 


4.2 Energy flow 

The following dehnitions are from [20) . We include them here for the sake of completeness. The 
analogous quantity of the Poynting vector in scalar wave theory is the energy flux vector; see 
|20j . We recall that for a real monochromatic field 

U{x,t) = Re [u{x)e~’'‘^^] , 

the averaged value of the energy flux vector, taken over an interval which is long compared to 
the period of the oscillations, is given by 

F{x) = —iC [n(x)VM(x) — n(x)Vu(x)] , 

where C is a positive constant depending on the polarization mode. In the transverse electric 
case, C = oj/pirn while in the transverse magnetic case C = ujjEm- We now consider the outward 
flow of energy through the sphere dB^ of radius R and center the origin: 

W = / F{x) ■ v{x)da{x), 

JdBa 

where v{x) is the outward normal at x G 

As the total field can be written asU = + tt®, the flow can be decomposed into three parts: 

W = W* + W" + W', 
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where 


w* = 

-iC 

/ 

u*(x)Vu*(x) 

— u*(x)Vu*(x) 

• b{x) da{x), 



IdBu 




w® = 

-iC 

[ [ 

u^{x)\7u^{x) 

— u^{x)Vu^{x)] ■ v{x) da{x), 



IdBu 




w' = 

-iC 

f 

u*(x)Vu^(x) 

— u^(x)Vu*(x) 

— tt*(x)Vu®(x 



IdBu 





In the case where n* is a plane wave, we can see that W* = 0: 


W* = - iC 


= -iC 


'OBr l 


tt*(x)Vrt*(x) — rt*(x)V«*(x) da{x), 


^^ikmd-x _|_ ^ikmd-xj^ ^^-ikmd-x 


ISBr 

=2Ckmd- / v{x)da{x), 
JoBr 


v{x) da{x), 


= 0 . 


In a non absorbing medium with non absorbing scatterer, W is equal to zero because the 
electromagnetic energy would be conserved by the scattering process. However, if the scatterer 
is an absorbing body, the conservation of energy gives the rate of absorption as 

= -W. 


Therefore, we have 


+ W® = -W. 

Here, W' is called the extinction rate. It is the rate at which the energy is removed by the 
scatterer from the illuminating plane wave, and it is the sum of the rate of absorption and the 
rate at which energy is scattered. 


4.3 Extinction, absorption, and scattering cross-sections and the optical the¬ 
orem 


Denote by C/* the quantity U''{x) = |u*(x)Vu*(x) — u*(x)Vu*(x) 
illumination, U'^(x) is independent of x and is given by 17* = 2krr 


In the case of a plane wave 


Definition 3. The scattering eross-section , the absorption cross-seetion and the extinc¬ 
tion cross-section are defined by 




jji’ 






qext 


-w' 

17* ■ 


Note that these quantities are independent of x. 
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Theorem 4.1 (Optical theorem). Ifu^{x) = where d is a unit direction, then 


km 

Q* = / \Aoo{x)\'^da{x) 


(4.3) 

(4.4) 


with Aao being the scattering amplitude defined by (4-2). 


Proof. The Sommerfeld radiation condition gives, for any x G dBn, 

Vu^(x) • ^{x) ~ ikmU^{x). 


(4.5) 


Hence, from (4.1) we get 


u^{x)Vu^{x) ■ v{x) — u^{x)Vu^{x) ■ v{x) 


2CK 


Aoo 


X 


\x\ 


which yields (4.4). We now compute the extinction rate. We have 

Vu''{x) ■ v{x) = ikmd ■ i'{x)e''^'^^'^. 


(4.6) 


Therefore, using 4.5 and 4.6 it follows that 


_ _ ^ikm{\x\-d-x) ^ikm.{\x\-d-x) 

u'^{x)Vu^{x) ■ v{x) — n^(x)Vti*(x) • v{x) =ikm -^—:—:- d ■ n — ik„ 


47r|a:| 
47 r|x| 


47r|x| 
{d ■ v{x) — 1). 


For X G dBn, we can write 


_ _ jU p—ikmRv{x)-{d—u{x)) 

u*(x)Vu®(x) • v{x) — u*(x)Vm*(x) • v{x) = — — --— - {d ■ v{x) — 1). 

47rii 

We now use Jones’ lemma (see, for instance, m) to write the following asymptotic expansion 
as i? —)• oo 

^ [ g(i/(x))e-*^-'^'^(*)dcj(x) ~ ^ (g{d)e-^^^^ - g{-d)P^'-^'] , 

RJdB„ km \ / 


where 


Hence, 


IdBn L 


g(iy(x)) = d ■ v{x) — 1. 


ri®(x)Vn^(x) — ti^(x)Vn*(x) • i/(x) ~ —Hoo(d) as 12 —)• oo. 
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Therefore 


W = -iC [A^id) - A^id)] = 2CQ [A^id)]. 


Since 


we get the result. 


\C\ u^{x)Vu\x)-u\x)Vui{x) =2\C\kr, 


□ 


4.4 The quasi-static limit 


We start by recalling the small volume expansion for the far-held. Let A be dehned by (2.19) 
and let 

(4.7) 


M{X,D):= [ {XId-IC*D)~^W]xda{x) 
JdD 


be the polarization tensor. The following asymptotic expansion holds. It can be proved by 
exactly the same arguments as those in [^. 

Propsition 4.2. Assume that D = 5B + z. As 5 goes to zero the scattered field can be written 
as follows: 


u%x) = -kf,{ - 1] \D\G{x,z,km)ufiz)-VzG{x,z,km) ■ M{X,D)Vufiz) 


+0 


dist(A, a{ICf))) 


(4.8) 


for X away from D. Here, dist(A, it(/C| 5 )) denotes miuj |A — Aj| with Xj being the eigenvalues of 


K 


D- 


Assume for simplicity that Sc = Sm- We explicitly compute the scattering ampli tude Aqo in 
(4.1). Take u^{x) = and assume again for simplicity that z = 0. Equation (4.8) yields, 

for \x\ > 


|3i| 


X X 


r-H ■ w J ■ ° Aist(A.a(/cB)) ^^ 

Since we are in the far-held region, we can write that, up to an error of order (^^/dist(A, a{ICf))), 

Jkm\x\ / ^ \ / 


u^{x) = —k^ 


A'k\x\ \|x 


— ■M{X,D)d] + 0 ^ . 


(4.9) 


In the next proposition we write the extinction and scattering cross-sections in terms of the 
polarization tensor. 

Propsition 4.3. The leading-order term (as 5 goes to zero) of the average over the orientation 
of the extinction cross-section of a randomly oriented nanoparticle is given by 


= ^Xs[TAM{X,D)] 


(4.10) 
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where Tr denotes the trace of a matrix. The leading-order term of the average over the orientation 
scattering cross-section of a randomly oriented nanoparticle is given by 


QL = ktn^\TTM{X,D)\^ 


(4.11) 


Proof. Remark from (4.9) that the scattering amplitude Aqo in the case of a plane wave illumi¬ 
nation is given by 




X 


^]=-kl^.M{X,D)d. 


(4.12) 


Using Theorem 4.1, we can see that for a given orientation 

= -km'^ [d ■ M(A, D)d] . 

Therefore, if we integrate over all illuminations we find that 


^ext _km^ 


L4S2 


d- M{X,D)dda{d) 


Since 'ijM(X, D) is symmetric, it can be written as QM{X, D) = P^N(X)P where P is unitary 
and N is diagonal and real. Then, by the change of variables d = P^x and using spherical 
coordinates, it follows that 


_km 


X ■ N{X)xda{x) 


= — [TriV(A)] = [TrM(A, D)] . 
3 3 


(4.13) 


Now, we compute the averaged scattering cross-section. Let iRM(X, D) = P*'N(X)P where P is 
unitary and N is diagonal and real. We have 


m —'“m II |x • M(A, D)d\ da{x) da{d) 

§2xS2 


D* =k 


§2xS2 


X ■ N{X)d da{x)da{d) -|- 


§2xS2 


x-N{X)d da{x)da{d) 


Then a straightforward computation in spherical coordinates gives 


Q^ = ^^^|TrM(A,Zl)|2, 


which completes the proof. 


□ 


From Theorem |4.1[ we obtain that the averaged absorption cross-section is given by 

|TrM(A,D)|^ 


kn 


Ql,= '^^[TTM{X,D)]-ktJf 


Therefore, under the condition (2.26), blows up at plasmonic resonances. 
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4.5 An upper bound for the averaged extinction cross-section 

The goal of this section is to derive an upper bound for the modulus of the averaged extinction 
cross-section of a randomly oriented nanoparticle. Recall that the entries D) of 

the polarization tensor M{\, D) are given by 


:= / xi{XI-JC*jj) ^[i^m]ix)da{x). 

JdD 

For a domain D in K,*j^ is compact and self-adjoint in T-L* (defined in Lemma 
d = 3 and in Lemma C.l for d = 2). Thus, we can write 


(4.14) 
for 


2.1 


{xid - icy) 


E' 

i=o 


A-A, 


with {Xj,(pj) being the eigenvalues and eigenvectors of in H* (see Lemma 2.1). Hence, the 
entries of the polarization tensor M can be decomposed as 


Mi,m{X,D) = Y, 


l,m 


i=i 


A - A. ’ 


(4.15) 


where •= ym,^j)n*yj^xi)_i i. Note that (j^m; x(^^))_i i = 0. So, considering the fact 

that Ao = 1/2, we have {um, V’o)n* = 0 and so, = 0. 

The following lemmas are useful for us. 


Lemma 4.1. We have 


Proof. For d = 3, we have 


j>i. 


{Pj:Xi)_i^i = Xj) \^Id-)Cy[<fj],xi'^_^ ^ 


-1 / dSpl^pj] 
1/2-AjV du 


,xi 


1 1 
■ 2 ’ 2 


dx, 


^^SD[pj]da - I (^AxiSdIpj] - xiASnlPjijdx 


JdD 

yi,Pj)n* 

1/2 - Xj ’ 

where we used the fact that is harmonic in D. The same result holds for d = 2 if we 

□ 


change Sd by Sd (see Appendix |C[). Since |Ajj < 1/2 for j > 1, we obtain the result. 
Lemma 4.2. Let 


Mi^yX,D) = Y, 


oo O') 


i=i 


A-A. 
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be the {l,m)-entry of the polarization tensor M associated with a domain D d Then, 
the following properties hold: 

(^) 


i=i 

(ii) 

j=i 1=1 

(Hi) 

CXD d f rj A\ d n 

^2 "!? = + ^ / 
j=i 1=1 1=1 

Proof. The proof can be found in Appendix □ 

Let X = X' + iX”. We have 

, , d |V|E?-ian 

|a(Tr(M(A.P)))| = E (y _ ('‘■16) 

For d = 2 the spectrum a{ICf))\{l/2} is symmetric. For d = 3 this is no longer true. 
Nevertheless, for our purposes, we can assume that (t(/C^)\{ 1/2} is symmetric by defining 
(i) 

afi' = 0 if Xj is not in the original spectrum. 

Without loss of generality we assume for ease of notation that Conditions!^ andhold. Then 
we define the bijection p : N"'' —)• N"'' such that A^q-) = —Xj and we can write 


1(TV(M(A,D)))| = X E 




+E 


\y\pipu)) 


^(A'-A,)2 + A"2 ^(A' + A,)2 + A"2 

|A"| “ (A'2 + A"2 + A|)(/3(j) + + 2X'Xj{P^^'l - 

~Y~ ^ 


i=i 


((A' - XjY + A"2) ((A' + XjY + A"2) 


where (dj = Eaff- 
z=i 

From Lemma Id.ll it follows that 

(A'2 + A"2 + A2)(/3(i) + /3 (pL))) + 2A'Aj(/3(^) - 
((A'-A,)2 + A"2)((A' + A,)2 + A"2) 


> 0 . 
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Moreover, 

(A'2 + A"2 + A2)(/50) + /3(^0))) + 2A'Aj(/3(^) - 
((A' - Aj)2 + A"2) ((A' + XjY + A"2) 


< 


(A'2 + A"2 + A2)(/3(j) + + 2A'Aj(/3(^') - 


A"2(4A'2 + A"2) 


Hence, 


|9(Tr(M(A,I))))|<^^ 


\//2 

^ MA'2 + A"2''- 


|A"| “ (A'2 + A"2 + A2)(/3(^') + + 2A'(Aj/3(^') + 


j=i 


A"2(4A'2 + A"2) 


+ o( 


A 


//2 


4A'2 + A"2^' 


Using Lemma 4.2 we obtain the following result. 


Theorem 4.2. Let M{X,D) be the polarization tensor associated with a C^’" domain D d 
with A = A' + iX” such that |A"| <C 1 and |A'| < 1/2. Then, 


KTr(M(A,Il)))| < 


+ 


d\y\\D\ 

A"2 + 4A'2 

1 

|A"|(A"2 + 4A'2) 
X!'^ 

^‘^^4A'2 +A"2) 


d>^‘^\D\ + ^^^\D\+Y, [ \VSDM^dx + 2X'^^^\D\ 
\ 1=1 dD 


The bound in the above theorem depends not only on the volume of the particle but also on 
its geometry. Nevertheless, we remark that, since |Aj| < 

i=i 1=1 

Hence, we can find a geometry independent, but not optimal, bound. 

Corollary 4.1. ITe have 

1 


KTr(M(A,Il)))| < 


|A"|(A"2 + 4A'2) 


^d\D\[X'^ + -) + 2A 


' '^A"2 + 4A'2^^^4A'2 + A"2^ 
(4.17) 
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4.5.1 Bound for ellipses 

If D is an ellipse whose semi-axes are on the xi- and X 2 - axes and of length a and b, respectively, 
then its polarization tensor takes the form [7] 


\D\ 


M{X,D) = 


\ _ 1 a—b 

^ 2 a+& 


0 

\D\ 


\ 


A + / 


On the other hand, it is known that in Ti* (dD) |27) 


(4.18) 




j = l,2. 


Then, from (4.15), we also have 


/ 


M{X,D) = 


E 


a 


U) 

1,1 


E 


a 


ii) 

1,2 


/=iA-E^V j=iA-E^') 


E 


2 a+b J 

U) 


a 


1,2 


E 


2 a+b J 

(i) 


a, 


2,2 




Let Ai = - 


1 a — 


2 a + b 


and V(Aj) = {i G N such that JC*^[^pi] = Xjipi}. It is clear now that 


E “S= E 


(i) 

« 2,2 = 


ieV(Ai) 


ieV(-Ai) 


for j >2 and 


for j > 1. 


1^1’ E "S = E “2^2 = 0 

ieV(Aj) ieV(-AT 


(4.19) 


E “L2 = 0 

ieV(Aj) 


In view of (4.19), we have 


/3(i) 


+ 


jjipij)) 


(A'-Aj)2 + A"2 (A' + Aj)2 + A"2 

Hence, 


^ 4A^2^(j) ^ Y/2 


A"2(4A'2 + A"2) 


4A'2 + A"2^' 


|il(Ti:(M(A,Zl)))|<^^ 


|A"| ^ 4A'2/3L) + A"2(/3(j) + /3(i)) 


j=i 


A"2(4A'2 + A"2) 


+ o( 


A"2 


4A'2 +A"2^- 


Using Lemma 4.2 we obtain the following result. 
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Figure 1: Optimal bound for ellipses. 


Corollary 4.2. For any ellipse D of semi-axes of length a and b, we have 

|5|4A'2 2|A"||5| A"2 

+ ,„o . ...o + 0 (- 


MTr{M{X,D)))\ < 


|A"|(A"2 + 4A'2) A"2 + 4A'2 MA'2 + A"2 


)■ 


(4.20) 


Figure shows (4.20) and the average extinction of two ellipses of semi-axis a and b, where 
the ratio afb = 2 and a/6 = 4, respectively. 

We can see from (4.16), Lemma 4.1 and the first sum rule in Lemma 4.2 that for an arbitrary 
shape B, |9(Tr(M(A, i?)))| is a convex combination of for Xj G it(/C^)\{ 1/2}. Since 

ellipses put all the weight of this convex combination in ±Ai = we have for any ellipse 

D and any shape B such that \B\ = \D\, 


\^{Tr{M{X*,B)))\ < |a(Tr(M(A*,L>)))| 


with X* = ±l^,+iX". 


Thus, bound (4.20) applies for any arbitrary shape B in dimension two. This implies that, for 
a given material and a given desired resonance frequency a;*, the optimal shape for the extinction 
resonance (in the quasi-static limit) is an ellipse of semi-axis a and 6 such that A'(a;*) = 


4.5.2 Bound for ellipsoids 

Let D be an ellipsoid given by 

The following holds [7]. 


^2 2 2 

■^1 I -^2 , -^3 _ -I 

Pi P2 Pi 


(4.21) 


Lemma 4.3. Let D be the ellipsoid defined by (4.21). Then, for x G D, 


SdWi]{x) = sixi, / = 1,2,3, 
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where 


P1P2P3 f°° _1_ 

2 Jo {pf + s)^/{pI + s){pi + s){pI + s) 


Then we have 


3 „ 

\VSD[i 2 i]\‘^dx = {sl + sl + sl)\D\. 

1=1 


For a rotated ellipsoid D = TZD with TZ being a rotation matrix, we have M (A, D) = IZM (A, D)7i?- 
and so Tr(M(A, D)) = Tr(M(A, D)). Therefore, for any ellipsoid D of semi-axes of length pi,P 2 
and po the following result holds. 


Corollary 4.3. For any ellipsoid D of semi-axes of length pi,p 2 andpo, we have 

.(TV(M(A,»))) s 1 °' Ai r ± . 4m ^ e( 


|A"|(A"2 + 4A'2) 


A"2 + 4A'2 


4A'2 + A"2 ' 


where for j = 1, 2,3, 


Si = -- 


P1P2P3 


lo [p] + s)^/{p{ + s){p ^2 + s){pi + s) 


ds. 


(4.22) 


5 Link with the scattering coefficients 

Our aim in this section is to exhibit the mechanism underlying plasmonic resonances in terms 
of the scattering coefficients corresponding to the nanoparticle. The concept of scattering co¬ 
efficients was first introduced in [9]. It plays a key role in constructing cloaking structures. It 
was extended in m to the full Maxwell equations. The scattering coefficients are simply the 


Fourier coefficients of the scattering amplitude ^oo- In Theorem 5.1 we provide an asymptotic 


expansion of the scattering amplitude in terms of the scattering coefficients of order ±1. Our 
formula shows that, under physical conditions, the scattering coefficients of orders ±1 are the 
only scattering coefficients inducing the scattering cross-section enhancement. For simplicity we 
only consider here the two-dimensional case. 


5.1 The notion of scattering coefficients 


From Graf’s addition formula [7j and (2.2) the following asymptotic formula holds as |x| —oo 


u^{x) = {u-u^){x) = 

n(PL 


\x\ e 


JdD 


where x = (|x|, 9x) in polar coordinates, Hn^ is the Hankel function of the first kind and order 


n, Jn is the Bessel function of order n and is the solution to (2.4). 
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we have 


For u^[x) = e 


ikmd-x 


"U (x) — ^ ^ Q'm(^ )'An(fcm|^|)^ 


imOx 


m£.7j 


where am{u^) = e *™^2 By the superposition principle, we get 

Ip = '^ 
m^'L 


where ipm is solution to (2.4) replacing / by 


fim) _ pim) ^ 

/ic ^ 


with 


1 


Ft\x) = - 


fdm 


du 


We have 


where 


U^{x) = {n- U^){x) = -lYl Y. 

nGlj mEZ 


Wnm = / Jnikm\y\)e *"^2'V’m(2/)c(cr(y)- 

JdD 


The coefficients Wnm are called the scattering coefficients. 
Lemma 5.1. In the space 'H*{dD), as oj goes to zero, we have 

/(±i) = ufYH0{u^), 

/M = 0{uj'^), \m\ > 1, 


where 


(±1) y/^mj-hn / 1 


/r ^=T 


1 A 


—+ -Qld - ICh)S^^[\x\A^^^ 


Proof. Recall that Jo{x) = 1 + O(x^). By virtue of the fact that 

(i/rf- (x:J,-)*)(s‘')-‘|x(8o)l = o(i.,2), 

we arrive at the estimate for (see Appendix . Moreover, 


J±i{x) = A- + 0{x^) 


(5.1) 
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together with the fact that 


(^/d- = {\ld-lC*i,)S^^ + 0{ojHogu:) 

gives the expansion of in terms of u (see Appendix . 

Finally, Jm{x) = 0{x'^) immediately yields the desired estimate for /I™). 

From Theorem |C.1[ it is easy to see that 

V’m = 2^- _ -h Ad(w) {Pjc{uj)f). 


j&J 


Tj{uj) 


□ 


(5.2) 


Hence, from the definition of the scattering coefficients. 


Wnm = Yl 

Since 


Jn{km\x\)e 


—inOa 


1 1 
■ 2 ’ 2 


Tj{uj) 


— + Jnikm\y\)e ^'^^yO{u:)da{y). 

IdD 

(5.3) 


Jm{x') 


ex \ l™l 


as m —)■ oo, we have 


y^27r|m|) V2|m| 

(j\m\ 


\f(m)\ < 


\m 


\m\ ■ 


Using the Cauchy-Schwarz inequality and Lemma 5.1, we obtain the following result. 
Propsition 5.1. For |n|, \m\ > 0, we have 

Q(^lJ\n\ + \m\) (j\n\+\m\ 

|72|hl|777,||m| 

for a positive constant C independent of lo. 

5.2 The leading-order term in the expansion of the scattering amplitude 

In the following, we analyze the first-order scattering coefficients. 

Lemma 5.2. Assume that Conditions 1 and 2 hold. Then, 
sr^ O(w^) „. . 

16 J ^ 


^±1 = E 

3&J 


+ 0{uj^\oguj) 


2 y fd'm f-^c 


Tjiuj) 


+ 0(a;). 
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Pro of. The expression of ipo follows from (5.2) and Lemma 5.1 Changing So by S p in Theorem 


2.1 


gives i{-Id — IC*jy)Sj^ [|x|e*^^], j =. Using now Lemma 


the expression of V'ii- 

Recall that in two dimensions, 


H 


5.1 


m 


yields 

□ 


UM = ^ ^ " -^)Xj + 0{u;Hoguj), 

where Xj is an eigenvalue of /C^ and Aq = 1 / 2 . Recall also that for 0 G J we need tj —)• 0 and so 
Hm —> oo, which is a limiting case that we can ignore. In practice, Pj{u})[cpo{uj)] = 0. We also 
have {ipj,x{dD))_i i = 0 for j 7 ^ 0. 

It follows then from the above lemmas and the expression (5.3) of the scattering coefficients that 

ur O log Ul) 

Woo — 2^ -UTTo- \-0{uj), 


j&J 


VLoii = ^ 




Tj{u) 

0{uj^ logo;) 

u(w) 


+ 0{uj), 




jeJ 


Tj{uj) 


Note that W±i±i has a special structure. Indeed, from Lemma|5.2| and equation (|5.3|), we have 
± ± w- 


3&J 

= E 

3&J 

_ 

4 


+ 0{u3^\ogU3) 




+ 0{U3 ), 


+ 0 (a;Mogu;) 

Tj{u3) 

E - - - - - + 0 ( 1 ) 


+ 0{u}‘^), 


k3&J 


A — Aj + 0{u!^ log w) 


where A is defined by (2.19). Now, assume that minjgj |Tj(a;)| S> w^logw. Then, 

o » o ^ 


= f (E 

\i6J 


A - A,- 


+ 0 ( 1 ) 


(5.4) 


Define the contracted polarization tensors by 

1V±,±(A,D):= [ |a;|e^*+(A/-/C| 5 )-^[e^*+](x)dfT(x). 
JdD 
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It is clear that 


N+,+ {\D) 
iV+,_(A,D) 
N-,+ {\D) 
N-,-{X,D) 


Mi,i(A,L») - M2,2{X,D)+i2M^4X,D), 
Mi,i{X,D) + M2,2{X,D), 

Mi,i(A,L») + M2,2(A,D), 

Mi,i(A, D) - M2,2(A, D) - i2Mi,2(A, D), 


where Mi^rn{X,D) is the (/,m)-entry of the polarization tensor given by ( |4.7[ ) 
Finally, considering the above we can state the following result. 


Theorem 5.1. Let Aoo be the scattering amplitude in the far-field defined in (fi2, 
incoming plane wave u*(x) = Assume Conditions 1 and 2 and 


min 

ieJ 


\Tj{oj)\ » logW. 


Then, A^^ admits the following asymptotic expansion 

X \ X T' 




x\ 


\x\ 


Wid + O{oj^), 


where 


IFi = 


^ f W-n + lFi_i - 21^1,1 i{Wi-i - W-ii) 

-iy_ii) -VF-ii - VFi_i - 2IFii 


Here, Wnm are the scattering coefficients defined by (5.1). 


Proof. From (4.12), we have 


X \ , 2 ^ ^ ^ 


— M{X,D)d. 


Since ICfi, is compact and self-adjoint in H*, we have 


Ni,±{X,D) = 

1=1 


A - A,- 


E--+ 0(1). 


16 J 


A - A.- 


We have then from (5.4) that 


-^iV+,+ (A,D) = W_ii + 0(a;2), 

-^iV+,_(A,41) = -Wn + 0(a;2), 
-^iV_,+ (A,Zl) = -Wn + Oicj^), 

-^N_^_{X,D) = Wi_i + 0 (a; 2 ). 


for the 
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In view of 


Mil = J(iV+,+ + iV_,_ + 2Ar+,_), 

M 22 = ^ (-A^+,+- A^-,-+ 2iV+,_), 

Mi 2 = ^(A^+,+ -iV_,_), 

we get the result. □ 

6 Super-resolution (super-focusing) by using plasmonic parti¬ 
cles 

It is known that the resolution limit (or the diffraction limit) in a general inhomogeneous space 
is determined by the imaginary part of the Green function in the associated space m- By 
modifying the homogeneous spaces with subwavelength resonators, we can introduce propagating 
subwavelength resonance modes to the space which encode subwavelength information in a 
neighborhood of the space embedded by the subwavelenghth resonators, thus yield a Green’s 
function whose imaginary part exhibits subwavelength peaks and therefore break the resolution 
limit (or diffraction limit) in the homogeneous space. The principle has been mathematically 
demonstrated in m- Here, using the fact that plasmonic particles are ideal subwavelength 
resonators, we consider the possibility of super-resolution (super-focusing) by using a system of 
identical plasmonic particles. The results in this section can be viewed as a consequence of the 
results in Section |3l 

6.1 Asymptotic expansion of the scattered field 

In order to illustrate the superfocusing phenomenon, we set 

^ik^\x-xo\ 

u\x) = G{x,Xo,km) = -y-u-T- 

47 r|x — xq\ 

Lemma 6.1. In the space H*{dD), as ui goes to zero, we have 

f = fo + 0 ( 1 x 62 ) -j- 0 ( 62 ), 


where fo = (/o,i, • • •, /o,l)^ with 
1 


fo,i = - 


A7:\zi - xoP 


(—{zi - Xo) ■ Tz{x) + —{\ld - K:\)Sj^][{zi - Xo) • (x - zi)] \ = 0(6'^ 
\hm he ^ / 


Proof. The proof is similar to that of Lemma 2.6 Recall that 


1 /I 


fi = Fi,2 + (/C^=/)(4‘;)-MTu]- 

fic ^ 
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We can show that 

1 1 5 3 

Fi ,2 = - ^- ^{Zi - xo) • i^ix) + 0{6-2 ) + 0{io5-2 ) in n*{dDi). 

nmol' 4:Trnm\zi - Xo\-^ 

Besides, 

ikm\zi-Xo\ 1 5 3 

u''{x)\aDi = -;rn- rXidDi) + — - ^{zi - xq) ■ {x - Zi) + 0(62) + 0{u)62) in n{dDi). 

4:Tr\zi — xo\ 4:Tr\zi — xo\ 

Using the identity — IC’^^)S'^^[x{dDi)] = 0, we obtain that 

This completes the proof of the lemma. □ 

We now derive an asymptotic expansion of the scattered held in an intermediate regime which 
is neither too close to the plasmonic particles nor too far away. More precisely, we consider the 
following domain 

Daa; = G min \x — zA ^6, max \x — zA <C 
’ i<i<L i<i<L k-’ 

Lemma 6.2. Let 'ipi G %*{dDi) and let v(x) = Lhen we have for x G 

v(x) = G(x,zi,k)(- — - — r-ik)T^ —W'/ 2/V’K2/)'^c^(y) + C>((5^)||V'/||^*(aD0 

^\x-zA '\x-zA JdDo 

+G{x,zi,k) / 'ipi{y)da{y). 

JdDo 

Moreover, the following estimates hold 


v{x) = 0(62) 

L 

v(x) = 0(62) 

L 


Proof. We only consider the case when 1 = 0 . The other case follows similarly or by coordinate 
translation. We have 


v{x) = S’l)['il>]{x) = [ G{x,y,k)'ilj{y)da{y) = - [ - -i){y)da{y) 

JdDo JdDo 47r|x - y| 

G(a:,y,t)=G(x,0,t)+ + E E 


Since 


|a=l| 




m>2 \a=m\ 


dy^ 
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and 


dG{x,0,k) _ 0 

dy» 47r|xnx| > \x\ \^y 

we obtain the required identity for the case I = 0. The estimate follows from the fact that 

Wy'^WnidDo) = o{6 2 ). 

This completes the proof of the lemma. 

Denote by 


□ 


Sj^i{x,k) = G{x,zi,k) 


X- Zl 

X-Zl\^ JqD^ 


yipj{y)da{y), 


Si{x,k) = G{x,zi,k) / (po{y)da{y), 

JdDo 

Hj,i{xo) = ^ - ^ 0 ) • ^^)H*{dDoy 

It is clear that the following size estimates hold 

Sj^i{x,k) = 0{6^), Si{x,k) = 0{6^), Hj^i{xo) = 0{6^) forj / 0, Ho,i{xo) = 0. 

Theorem 6.1. Under Conditions^^^ an the Green function r(x, xq, km) in the presence 
of L plasmonic particles has the following representation in the quasi-static regime: for x G 
^5,km ; 

T{x,XQ,km) = G{x,xo,km) 

_l_ ^j,pi^o)^j,i,p^j,hg^j,gi^^ ^rti) + ^ ^ ^3 

jeJ 1=1 + 'DV + 0 {S‘^) + 0 (< 52 w 2 ) 

Proof. With n*(x) = G(x,xo, km), we have 

'^ = '^ Y1 + 0{6^, 

jej i<i<L i<i<L 

where 


0 -j,l — if, Pj,l)H*{dD) — {fo, Pj,l)n*{dD) + 0(1062) + 0(62), 

= (- )Xj^i^pHj^p(xo) + 0(0:62) + 0(62), 

f^c l^m 

Oo,/ = if,P0,l)'H*{dD)=O(62). 
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By Lemma 6.2 


l<p<L 1<P<L 

= ^ fern) + 0(5^) + 0(w5i). 


l<p<L 

On the other hand, for j = 0, we have 

?k. 


•Sir[v’o,i]{x) = 0 { 62 ), 

ro,iM = ro + 0(0 + 0(^V) = 0(l). 


Therefore, we can deduce that 


u = 


SD^mx) = Yl E E ao,i5^’"Kz]+0(<5"), 

jeJ i<z<i i<«<i 

Z; 1 - - 

y'y' - )^jA^o)^j,l,p^j,l,qSj,qix,km) + oM^) + 0{6- 

^ k'C fj^m 

+0{5^), 

Hj,p{xo)Xj^i^pXj^i^qSj^q{x, km) + 0{uj6^) + C>((5'^) ^ q/-j3\ 
jej I=i + {-j^- ]^) Tj,i + 0(<5^) + 0{6‘^u‘^) 


6.2 Asymptotic expansion of the imaginary part of the Green function 


□ 


As a consequence of Theorem 6.1, we obtain the following result on the imaginary part of the 
Green function. 


Theorem 6.2. Assume the same conditions as in Theorem 6.1. Under the additional assump¬ 
tion that 

A-Aj + (—> 0 { 6 ^) + 0 { 5 ^uj^), 

he Tm 


for each I and j G J, we have 

^T{x,xo,km) = ^G{x,xo,km) + 0(6^) + 

L 

EE 3f? (ifj,p(xo)A,y,pA,-,(x, 0) + OM^) + 0(5^)^ 

1 \ 


jEJ 1—1 


+ Si. ’ 
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where x,xo € 


Note that 5R \^Hj^p{xo)Xj^i^pXj^i^qSj^q{x,0)j = 0{6^). Under the conditions in Theorem 
if we have additionally that 


6.2 


+ {h ~ iL) 


-1 




= 0(;p) 


for some plasmonic frequency oj, then the term in the expansion of 3‘r(x, xQ) ^m) which is due 
to resonance has size one and exhibits subwavelength peak with width of order one. This breaks 
the diffraction limit 1/km in the free space. We also note that the term QG{x,xo,km) has size 
0(w). Thus, we can conclude that super-resolution (super-focusing) can indeed be achieved by 
using a system of plasmonic particles. 


7 Concluding remarks 

In this paper, based on perturbation arguments, we studied the scattering by plasmonic nanopar¬ 
ticles when the frequency is close to a resonant frequency. We have shown that plasmon resonant 
nanoparticles provide a possible way not only of super-resolved imaging but also of scattering 
and absorption enhancements. 

We have derived the shift and broadening of the plasmon resonance with changes in size. 
We have also consider the case of multiple nanoparticles under the weak interaction assumption. 
The localization algorithms developed in [ZIEIIIS] can be extended to the problem of imaging 
plasmonic nanoparticles. We have precisely quantified the scattering and absorption cross- 
section enhancements and gave optimal bounds on the enhancement factors. We have also 
linked the plasmonic resonances to the scattering coefficients and showed that the leading-order 
term of the scattering amplitude can be expressed in terms of the ±-one order of the scattering 
coefficients. 

The generalization to the full Maxwell equations of the methods and results of the paper is 
under consideration and will be reported elsewhere. Another challenging problem is to optimize 
the super-focusing phenomenon in terms of the organization of the nanoparticles. This will be 
also the subject of a forthcoming publication. 

A Asymptotic expansion of the integral operators: single par¬ 
ticle 

In this section, we derive asymptotic expansions with respect to k of some boundary integral 
operators defined on the boundary of a bounded and simply connected smooth domain D in 
dimension three whose size is of order one. 

We first consider the single layer potential 

S%W\{x)= ( G{x,y,k)'ip{y)da{y), x € dD, 

JdD 
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where 


G{x,y,k) = 


^ik\x-y\ 

47r|x — y\ 


is the Green function of Helmholtz equation in subject to the Sommerfeld radiation condition. 
Note that 


(^( ^ iik\x-y\y _ 1 ik ^ {ik\x - y\y ^ 

j\4:TT\x — y\ 4 :Tt\x — y\ j! 


We get 

OO 

i=i 


where 

= - i^{y)dcr{y). 

In particular, we have 


(A.l) 


•SnAi’Kx) = f ip{y)da{y), 

47r JdD 

SD, 2 [ij]{x) = -— / \x - y\il!{y)da{y). 

47r JdD 


(A.2) 

(A.3) 


Lemma A.l. \\<SD,j\\c{{'H*{dD),'H{dD)) uniformly bounded with respect to j. Moreover, the 
series in (A.l) is convergent in C{(H*{dD),'H{dD)). 


Proof. It is clear that 


\\'SD,j\\c{L^idD),H^{aD)) < C, 


where C is independent of j. On the other hand, a similar estimate also holds for the operator 
Sf) J. It follows that 

\\'SD,j\\c{H-paD),L2{dD)) < C. 

Thus, we can conclude that ||5n Jl i , i , is uniformly bounded by using interpo- 

lation theory. By the equivalence of norms in the H~'i{dD) and H^^dD), the lemma follows 
immediately. □ 

Note that 5£) is invertible in dimension three, so is for small k. By formally writing 

(5^)-^ = 5^1 + kBo,! + k^BD,2 + ..., (A.4) 

and using the identity = Id, we can derive that 

Sd,i =—•Sf)^SD,i'Sf)^, Bd, 2 =—Sf)^SD,2'Sf)^ + Sf)^SD,iSf)^SD,i'Sf)^. (A.5) 

We can also derive other lower-order terms Bdj- 


Lemma A.2. The series in {A. 4 ) converges in C{TL{dD),TL*{dD)) for sufficiently small k. 
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Proof. The proof can be deduced from the identity 


i=i 


□ 


We now consider the expansion for the boundary integral operator (/C^)*. We have 


(/C^)* =fCh + kfCD,i + k^fCD,2 + ..., 


(A.6) 


where 

In particular, we have 


^D,i = 0, 


JCd, 2 [^P]{x) = ]- [ 
4vr JdD 


{x -y) ■ lyjx) 
\x - y\ 


f’{y)d(T{y). 


(A.7) 


Lemma A.3. The norm \\K,D,j\\c(H*{dD),H*{dD)) uniformly bounded for j > 1. Moreover, the 
series in ( A.6) is convergent in C(fH*{dD),TL*{dD)). 


B Asymptotic expansion of the integral operators: multiple 
particles 

In this section, we consider the three-dimensional case. We assume that the particles have size 
of order 5 which is a small number and the distance between them is of order one. We write 
Dj = Zj + 5D, j = 1,2,..., M, where D has size one and is centered at the origin. Our goal 
is to derive estimates for various boundary integral operators considered in the paper that are 
defined on small particles in terms of their size. For this purpose, we denote by Dq = 5D. For 
each function / defined on ODq, we define a corresponding function on D by 

V{f){x) = f{dx). 


We first state some useful results. 

Lemma B.l. The following scaling properties hold: 

(i) h{f)\\L2(dD) = d~^\\f\\L^dDo); 

(ii) \\vif)\\u(dD) = d~M\f\\H{dDo); 

(iii) Mf)\\u*(dD) = d~h\f\\H*{dDo)- 
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Proof. The proof of (i) is straightforward and we only need to prove (ii) and (hi). To prove (iii) 
we have 


H*(dDo) 


I I 

JdDn J 81 


f{x)f{y) 


= 5 


dDo JdDo 4vr|x - y\ 

3 / / V{f)ix)r]{f){y) 


IdbJdD 47r|x-y| 


da{x)da{y) 

daix)da{x) 


whence (iii) follows. To prove (ii), recall that 

WfWuidDo) = II‘5£,q/||-h*(9Z)o)' 

Let u = Sfy^[f]. Then f = SdoIu]- We can show that 

Vif) = 

As a result, we have 

\\dif)\\n{ab) = ^W-SbidiuMuidb) = ^Mu)\\u-{ab) = ^~^\MH*{aDo) = ^~^\\f\\H{aDo), 

which proves (ii). □ 

Lemma B.2. Let X and Y be bounded and simply connected smooth domains in Assume 
0 G X,Y and X = 6X, Y = 5Y. Let TZ and TZ be two boundary integral operators from V'{dY) 
to V'{dX) and V'{dY) to V'{dX), respectively. Here, V denotes the Schwartz space. Assume 
that both operators have the same Schwartz kernel R with the following homogeneous scaling 
property 

R{6x,6y) = 5^R{x,y). 

Then, 

\\'T^\\c('H*(aY),-H*(ax)) = \\^\\c{'H-(ay),'n*{ax))-’ 

\\'^\\c('H*(aY),-H(ax)) = \\^\\c{'H-(aY),'H.{ax)y 

Proof. The result follows from Lemma |B.1| and the following identity 

oTZoy. 

□ 


We first consider the operators S^, and {Kff),)*. The following asymptotic expansions hold. 
Lemma B.3. (i) Regarded as operators from H*{dDj) into 'H{dDj), we have 

S% = Sd, + kSD,,i + k^SD„2 + 
where Sdj = 0(1) and SDj,m = 0(5™); 
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(a) Regarded as operators from T-L^dDj) into 'H*{dDj), we have 

{S’hy = s^] + kBo^i + k^BD„2 + 0{k^5^), 

where = 0(1) and Boj^m = 0(6"^); 

(Hi) Regarded as operators from 'H*{dDj) into 'H*{dDj), we have 

{K%r = lC*j,^+k^0{5^), 


where IC*^, = 0(1). 

Proof. The proof immediately follows from Lemmas |B.2[ |A.1[ and A.3 
We now consider the operator . By definition, 


Sd,,DiW\{x)= I G{x,y,k)'tlj{y)da{y), x e dDi. 

" JdDi 


Using the expansion 

where 

we can derive that 

where 

We can further write 

where is defined by 


oo 


G{x,y,k) = ^ k'^Qm{x,y), 

m=0 

^m\x_y\m-l 


Qm{x, y) — 


dvr 


m>0 


•Sj,i,m[il^]{x) = / Qm{x,y)i>{y)da{y). 

JdDi 


Sj,l,m — ‘S'jY 


Lm.ri’i 


n>0 


1 5l«l+l/3| 

2 ^ PPhd7P:d^BQrn{zuZj){x-Zir{y-Z^ 


dDj I I dx<^dyh 

J \a\ + \y\=n 


□ 


?)^V’(y) da{y). 
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In particular, we have 




— ~47r|2;- - X(^-^i))-H'-l/2(aDj),_H'i/2(ai5^)X(A), 

= }^L 1^1 13 X{^Dl))H-l/2^QD^)^HU2^^9D.^ + ((y - ^j)“, 

|a|=l 



i«i+i/3i=2 

'I 

= “^(V’l x{dDj))jj-i/2(^Q]:).'jjji/2(^Qjy.'jx{Di) 


Sj^l^l[0\{X] 


Sj ^1^2,oi'^Kx) = 

The following estimate holds 


\zi - Zjli'ip^xidDj)) H-^/^{dDj),H^/^{dDj)X{Di). 


ine toliowmg estimate nolds. 

Lemma B.4. ITe have \\Sj^i^jn,n \\c{H*(80)^90)) < 

Proof. After a translation of coordinates, the stated estimate immediately follows from Lemma 

[Rl □ 

Similarly, for the operator /C^ defined in the following way 


l<^D,,Di['<P]ix)= [ ^^^wY^fj{y)dcr{y), X 

JdDi Ol2[X) 


edDi, 


we have 


y'k 

^D,,n 


m>0 n>0 


,m,n? 


where 

ir \ n/ \ — f ^ d’^Kmizi-, 

- J Y Q^hQy, 

J |o| + |/3|=n 




with 


^ ^ i™'(m — l)|x — yp ^ 

A,^(x,y) = . 


In particular, we have 

/Cj'Z,o,o[V'](^) = ^^.|3 • ^(^)(V^^X(g^i))H-i/20D,).gV20DT 

— (V') (y — • ^i^)) H-^/' 2 {dDj),H^/^{dDj) 

+ izi — Zj) ■ V'{x)['ll^, x{dDj)) f^-l/ 2 (^g£,.^ jjl/ 2 (^QD^'^ , 
J<^j,l,i,m['ip] = 0 for all m. 

We have \\ICj^l^rn,n\\c{n* (aDj),H*{ODi)) ~ 


(B.l) 

(B.2) 


Lemma B.5. 
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Proof. Note that 


^j,l,m,n [^] (®) 



1 d^KrajzuZj) 

\j3\ dx^dy^ 

1 d’^KmjzuZj) 

a!/3! dx^dy’^ 

1 d^Km{zi,Zj) 

a!/3! dxPdy°‘ 


{x 


{x 


zif{y - ZjY{x - zi) ■ v{x)ii{y)da{y), 

- ziY{y - ZjY{y - Zj) ■ v{x)ii{y)da{y), 

- ziY{y - ZjY{zi - Zj) ■ v{x)'ip{y)da{y). 


After a translation of coordinates, we can apply Lemma [B. 2 to each one of the three terms above 
to conclude that /Cjy,m,n = This completes the proof of the lemma. □ 


To summarize, we have proven the following results. 


Lemma B.6. (i) Regarded as an operator from 'H*{dDj) into T-L{dDi) we have, 

‘^Dj,Di ~ ‘^j,1,0,0 + 0,2 + + hf^Sjj^2,Q + O(h^) + 0{k‘^5‘^). 

Moreover, 

Sj,l,^,n = 0{5^+^). 


(a) Regarded as an operator from PL* [dDj) into Pi*{dDi), we have 

JC%^n,=f^j,i,o,o + 0{k^6^). 

Moreover, 

1,0,0 = 0{5^). 


C Adaptation of results to the two-dimensional case 


In this section we adapt the layer potential formulation to plasmonic resonances in two dimen¬ 
sions. We only consider the single particle case. For the multiple particle case, a similar analysis 
holds. 

Recall that in the single-layer potential Sd '■ H~^/‘^{dD) —)• is not, in general, 

invertible nor injective. Hence, —(n,5D[n])_i i does not dehne an inner product and the sym- 


metrization technique described in Lemma 2.1 is no longer valid. 

To overcome this difficulty, a substitute oi Sd can be introduced as in 13] by 


SnVf] = 


SdW\ if = 0, 

1 if V' = Po, 


(C.l) 


where ipQ is the unique (in the case of a single particle) eigenfunction of IC*^ associated with 

eigenvalue 1/2 such that ((/^o, x(9T>))_ i i = 1. Note that, from the jump relations of the layer 

2 ’ 2 

potentials, is constant. 

The operator Sd ■ H~^/'^{dD) —)• is invertible. Moreover, the following Calderon 
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identity holds KLdSd = SdK,*j^. With this, define 


(u,v)n* = -{u,Sd[v])_ii. 


Thanks to the invertibility and positivity of —Sd, this defines an inner product for which is 
self-adjoint and Ti* is equivalent to Then, if D is C^’“, we have the following result. 


Lemma C.l. Let D be a hounded simply connected domain ofM? and let Sd be the operator 


defined in C.l. Then, 


(i) The operator K.*^ is eompact self-adjoint in the Hilbert space 'H*{dD) equipped with the 
inner product defined by 

{u,v)n* = -{u,Sd[v])_ii (C. 2 ) 

with (•, •)_! 1 being the duality pairing between and H^/^idD), which is equiv- 

2 ’ 2 

alent to the original one; 


(a) Let {\j,ipj), j = 0,1,2,..., be the eigenvalue and normalized eigenfunction pair of JCf) 
with Ao = 5 . Then, Xj G 5 ] o,nd Xj —>■ 0 as j —)• 00 ; 

(in) T-L*{dD) = T-lQ{dD) Q {/aipQ, p, G C}, where 'HQ{dD) is the zero mean subspace of'H*{dD); 

(iv) The following representation formula holds: for any tf G , 


= X] ® Tj- 

j=0 


Remark C.l. Note that Sj^^[x{dD)] = ipo and {—^Id-\-}Cf)) = {—^Ld-\-JC)^)V'^*, where V-n* 
is the orthogonal projection onto ^{^{dD). In particular, we have {—^Id-\- fCf))Sfi)^[x{dD)] = 0. 

Let us now consider the single-layer potential for the Helmholtz equation in 

•SDbP]{x)=[ G{x,y,k)-fi{y)da{y), x e dD, 

Jan 


where G{x,y,k) = —^HQ^\k\x 
We have 


y|) and 


is the Hankel function of first kind and order 0. 


i 

4 


Hl^^\k\x 


where 


y\) 


1 


log |x 


y\ + Xk + '^{bj log k\x-y\ -ICj){k\x 
i=i 




Tk = 7i-(loS^ + 7 - Iog2) 

ZTT 



27r 22i(j!)2’ 


-bj 


^7 - log 2 


iTT 

~2 


E 

n=l 


1 

n 
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and 7 is the Euler constant. Thus, we get 


5^ = 5^ + E ^DJ + E ’ 

i=i i=i 


(C.3) 


where 


shmx) 

= SD[tp]{x)+Tk / ['tp]dcr, 


JdD 

‘527[V’](a;) 

= / bj\x - y\‘^^'ip{y)da{y), 

JdD 


= / \x-y\‘^^ibjlog\x-y\+Cj)'il>{y)da{y) 

JdD 


Lemma C.2. The norms WS^^WcW*00) ^00)) II‘5_d jllr:('H*(5r>),'H(9r>)) are uniformly bounded 
with respect to j. Moreover, the series in ( ’(7.5] j is convergent in £{'H*{dD),'H{dD)). 

Proof. The proof is similar to that of Lemma |A.1[ □ 

Observe that 

^Sd - Sd) [V’] = (•Sd -Sd) ['Pn*^^ + if^,To)n*To] = ii’,To)n* - x{dD)). 


( 2 ) 


Then it follows that 

•SdW = + i-ip, To)n* (‘5ri[</?o] - x{dD)) + Tk [ V’o + (V', (po)n*Todo- = 5d['0] + TkbP], 

JdD 

where 

'^ki'ip] = To)h* (5d[v2o] - x{dD) + Tk) ■ (C.4) 

Therefore, we arrive at the following result. 

Lemma C.3. For k small enough : %*{dD) — ^{{dD) is invertible. 

Proof. Tk is clearly a compact operator. Since Sd is invertible, the invertibility of S^ is equiv¬ 
alent to that of = Id + TkSf^^ . By the Fredholm alternative we only need to prove the 

injectivity of Id -|- 

Since V n G G C, for (^Id + T^Sf^^^ [n] = 0, we need v = Solot^n] = a G C. 

We have 


[id + TkSf,^) 




= a(5D[v?o] + T-fc) = 0 


iff 5 d[(/?o] = -Tk or a = 0. 


Since we can always find a small enough k such that 7^ —Tk, we need a = 0. This yields 

the stated result. □ 

Lemma C.4. For k small enough, the operator : 'H*{dD) —)• 'H{dD) is invertible. 
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Proof. The operator is a compact operator. Because is invertible for k small enough, 

by the Fredholm alternative only the injectivity of is necessary. From the uniqueness of a 
solution to the Helmholtz equation we get the result. 


□ 


We can write (C.3) as 


+ Gk^ 


where Gk = k'^ log kS^j^\ + + 0{k^ log k). From the two lemmas above we get the identity 


-1 


It is clear that \\{<S^)~^\\c{'H{dD),'H*{dD)) is bounded in k. Thus, for k small enough, we can 
formally write 

- {S^^r^GkiSl^r^ + OikHogH). 


We have the identity 




ck 


-1 


— 1 


D ■ 


Here, 

Then, 

and therefore. 


Kk = Id + (•,<Fo)w*(5d[</?o] -x{dD) +Tk)<.po. 


K = Id - {■,(Po)h* 


-x{dD) + Tk 




0%) ^ — {Sj^^[-]^ipQ)'kL*ipQ + 


5d[v^o] + Tk 


5ri[v7o] +Tfc 


Po- 


Finally, we get 

= CD+Uk-k"^ log kCDS^^\CD - k"^ (CdS^^\Cd - log k{UkS^^\CD + 


+0(A:^log ^ k) 


with Cd = and Uk = 


{^D [•])£o)'H* 


(pQ. We note that Uk = 0(log ^ k). 


■■o ^ SDVPo] + Tk 

We now consider the expansion for the boundary integral operator (/C^)*. We have 

OO OO 

(/C^)* = /Cl, + £ log k) + £ k^^lC% 
i=i j=i 


(C.5) 
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where 


'4){y)da{y), 

JdD av{x} 




^ dv{x 

r d {\x-y\‘^^{bj log\x-y\ + cj)) 

IdD ^{x) 


i^iy)dcr{y). 


Lemma C.5. The norms \\}C^nJ\c(n*(aD),n*(p)) and \\}C^D,j\\c{n*{dD),n*idD)) are uniformly 
bounded for j > 1. Moreover, the series in (C.5) is convergent in C{TL*{dD),TL*{dD)). 


Proof. The proof is similar to that of Lemma A.l 


□ 


Recall (2.5) and (|2.6|), we can show that the following result holds. 


Lemma C.6. Regarding Ad{oj) as an operator from T-L*{dD) to T-L*{dD), we have 

Aoiio) = Ad,o + w^(loga;)^D,i + O(w^), 

where 

-^D,o = 




Ad,1 = {Srald - ScVui) + —Qld - )C}))S {pmSmld - RcSc'Ph^o) ■ 


1 .1 


— 1 c(l) 


he 2 


Proof. We have 




+ hlk, - Up‘{logUj)£chcddDS^^\CD + 0{u?) 
= Sd + + Lo'^ {logUj)£mhm‘S^^i + O(w^). 


Also, CoAkrr, = V-hi{Sd) = 0, where Tfc^ is defined by ( |C.4[ ). Hence, 

(‘^ d ) = V'H* +lAk^SD Aldk^Tk^ + lJ^{\ogLLl)(£mhmT-DS\^'^l — £chcddDS\^'^lCDSD) + 0 {lJ^) 


— 'P'Hl AUk^So T logioCD<S\j^i[£mhmId — ScheP-w^) +0{uj^). 


From Remark C.l, it follows that 


{ Id-lC*j,)Uk^ = Q. 


Since \ld — (/C^)* = {\ld — /C^) — log + 0{k‘^), we get the desired result. 

Under Conditionsandthe perturbed eigenvalues and eigenfunctions of Ad{<jj) have the 
following form 


□ 


Tj{u}) = Tj + w^(log w)rj -4 + 0(a;^), 
(Pj{co) = (f j + ui'^ {log co)(pj^i + 0{oj'^), 


(C.6) 

(C.7) 
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where 


Tj.l — Rjj, 




Rji 




and 


Rjl = {•^D,l[^j],^l)H*- 


(C.8) 

(C.9) 


It is clear that Lemma 2.5 holds in the two-dimensional case. We also have the following 
asymptotic expansion for / in terms of ui. 

Lemma C.7. In the space R*{dD), as u goes to zero, we have 

f = ujfi + O(uj^), 

where 


fi = (-[d- izix)] + —{^rld - IC}j)Sj^^[d ■ (x - z)] 

^ h-m he ^ 


and z is the center of the domain D. 

Finally, the following result holds. 

Theorem C.l. Under Conditions'^ an the scattered field by a single plasmonic particle, 
= u — vI, has in the quasi-static limit the following representation: 


where 




^ikmU^'^'^'^[d-v{x),ipj)^,ipj + 0{uj'^\ogu:) 

V’ = 2^-^^ _^ -+ o{u) 


ieJ 


\ 0(a;2 logw) 


with A being defined by (2.19). 
Proof. We have 


; ^ ^-l^D : \f\ 

= 2 ^-- +Ad{uj) {Pje{u)f) 


j&J 


Tjiw) 


Y- {f 1, Pj)+ 0 {u:^{loguj)) 

1 . 1 - 7U - 1 9, ^ + 0(w). 


,6J Sii + 2^ - log 


IV) 


Since d- {x — z) is a harmonic function, changing So by Sd in Theorem 2.1 yields 
^/d-/CL)6ndd-(x-z)|,(p, , 

Then, the proof is complete. 


({\dd - lCj:,)Sj^\d ■ {x - z)],ipj'^^^ = -{d ■ iy{x), 


□ 
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Corollary C.l. Assume the same conditions as in Theorem 2.1. Then, under the additional 
condition 


min |Tj(a;)| ;:?> to , 
3&J 


we have 


^ +0(a;3logw) 

^ = 1 . —^ ./I — rw -+ 


ieJ 


A-Aj+w2loga;(^-^) Tj,i 


D Sum rules for the polarization tensor 

Let / be a holomorphic function defined in an open set 17 C C containing the spectrum of JCgj^ 

OO 

Then, we can write f{z) = for every z G U. 

j=0 


Definition 4. Let 


mh) :=Y.a,{JChy, 

j=0 


where := fC^ o fC*^ o .. o fC*jg . 

'-V-' 

j times 

Lemma D.l. ITe have 


f{ICh) = Y.f{Xj)i;Tj)n*Tj- 

i=i 


Proof. We have 


OO OO 


filCh) = 


Y^^iilChY = 

i=0 i=0 j=l 


OO / OO 


I (■) 

j=l \i=0 / 

OO 

j=i 


From Lemma D.l we can deduce that 

/. OO 

xif{lCjg)[v.m\{x)da{x) = '^ f {Xj)a\^l^. 


IdD 


i=i 


□ 


(D.l) 


Equation (D.l) yields the summation rules for the entries of the polarization tensor. 
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oo 

In order to prove that = ^i,m\D\, we take /(A) 

i=i 


1 in (D.l) to get 


oo „ 

m = / XlVm{x) da{x) = 5l^r. 

,=1 ’ 


II?I 


Next, we prove that 


OO d 


j=i «=i 


(i) _ 


£>l 


Taking /(A) = A in (D.l), we obtain 

OO d 


Yj Y “m ^Y (^) 

j=i i=i i=i 

[ xiJC*j^[i^i]{x) da{x) = [ xi (Yi{x) + 

Ian JdD 

\D 


2’"-"' ' du 

dSD[n\{x 


da{x), 




da{x). 


(D.2) 


Integrating by parts we arrive at 
dSoWiKx 


Ud'''' du 


{x)da{x) = / ei{x) ■'S/SD[i'i]{x)dx + / xiASD[j^i]{x)dx. 

Jd Jd 


Since the single-layer potential is harmonic on D, 
dSDWi]{x 


L""'” diy 


{x)da{x) = / ei{x) 
JD 


'dD 


Vx^{x,x')i'i{x')da{x') \ dx. 


Summing on i and using Va:r(x,x') = —Vx'^{x, x'), we get 

Y[ ^du{x) i^)d(j{x) = -J v{x) ■Vx''^{x,x')da{x'^ dx, 


1=1 


IdD 


= - 'DD[l]{x)dx, 
Jd 

= - IDL 


(D.3) 


where Pd is the double-layer potential. Hence, summing equation (D.2) for i = 1,...,d, we get 
the result. 

Finally, we show that 

OO d y-J A d n 

Y^^iY^U = ^“ 1^1 [ l^SoWifdx. 


j=l 1=1 


1=1 


ID 
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Taking /(A) = in (D.l) yields 


CO d d 

„(i) 


Y “m =Y da{a 

j=i 1=1 '-1 


i=i 

d 


Y / X^D[yi]ix)IC})[ui]{x)da{x) 


i=i 


IdD 


t,f + j: f_IColy,]^^\-d. 

l=l ddD 


IdD 




yi QSdWi] 


du 

d 


l=l 


'dD 


2 du 


dfx + Y i 


dSnWi] 


du 


da. 


h 


h 


From (D.3) it follows that 




Since x; is harmonic, we have xi = 'DD[yi]{x)\- — SD[i^i]ix) on dD, and thus, 




i=i 


IdD 


da{x), 




l=l 

d 


IdD 


da, 


D\ + Y \^SDh]fdx. 

i=i 


Replacing Ii and I 2 by their expressions gives the desired result. 
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